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Context. Apparent period variations detected in several eclipsing, close-compact binaries are frequently interpreted as 
being caused by circumbinary giant planets. This interpretation raises the question of the origin of the potential planets 
that must have either formed in the primordial circumbinary disk, together with the host binary star, and survived its 
^Lj ■ evolution into a close-compact binary or formed in a post-common-envelope circumbinary disk that remained bound to 

the post-common-envelope binary (PCEB). 

Aims. Here we combine current knowledge of planet formation and the statistics of giant planets around primordial 
and evolved binary stars with the theory of close-compact binary star evolution aiming to derive new constraints on 
possible formation scenarios. 

Methods. We compiled a comprehensive list of observed eclipsing PCEBs, estimated the fraction of systems showing 
apparent period variations, reconstructed the evolutionary history of the PCEBs, and performed binary population 
models of PCEBs to characterize their main sequence binary progenitors. We reviewed the currently available constraints 
C"| on the fraction of PCEB progenitors that host circumbinary giant planets. 

O [ Results. We find that the progenitors of PCEBs are very unlikely to be frequent hosts of giant planets (< 10 per cent), 

while the frequency of PCEBs with observed apparent period variations is very high (~ 90 per cent). 
Conclusions. The variations in eclipse timings measured in eclipsing PCEBs are probably not caused by first-generation 
,£j . planets that survived common-envelope evolution. The remaining options for explaining the observed period variations 

are second-generation planet formation or perhaps variations in the shape of a magnetically active secondary star. We 
suggest observational tests for both options. 
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ly-) , 1. Introduction Circumbinary exoplanets around binary stars composed of 

rf) „ , , , , n • , t w0 MS stars (MS+MS binaries) were predicted long ago, 

The discovery of extrasolar planets (exoplanets) m the but only recently have the first six d etections been possi ble 
nineties and the subs equent detection of m ore than 700 ggg^ 2011 Wel sh et al.ll2012llOrosz et alfeoillbl ). 
exoplanets until now (pchneider et al.l 120111 ) has taken us 
- much closer toward answering fundamental questions, such Amon S the Perhaps most unexpected potential planet 
as how planets form from circumstellar material and what detections are those derived from apparent period varia- 
the fate of planetary systems is when their host stars evolve tlons ( or ecll P se tlme variations), frequently measured m 
into giant stars and compact objects. Especially that exo- eclipsing, close-compact binary stars containing a white 
planets have been detected in rather unexpected exotic en- dwarf ( WD ) or hot subdwarf B (sdB) primary star and 
vironments might go further pave our way towards solving a low-mass star companion. The first of these substellar 
n3 ' these important questions third bodies, a brown dwarf (BD), was announced roughly 

Pulsar timings have led to detecting the first con- a decade a S° around the post-common-envelope binary 
firmed exoplanet orbiti ng the pu lsar PSRB1257+12 (PCEB) V471 Tau (|Guman fc Ribas| | 200H ). At present, de- 
ll Wolszczan fc Fraill[l992l IWolszczanl Il99l . and continued tectl0n of substellar and mostly planetary circumbinary 
monitoring of the system has reveale d the existence of a companions has been claimed for a dozen eclipsing, close- 
second and a third planet (lWolszczanlll994D . even the pos- compact binaries containing WD or hot sdB primary stars 
sibility of a fourth one (|Wolszczanl llflolT The first exo- (see | Beuermann et al. ||201fl for the most convincing ex- 
planet around a solar-type star, a Jupiter-mass companion am P le )- Whlle the Parameters derived for these potential 
to the star 51 Pegasi, was id entified and confir med thanks circumbinary giant planets have to be considered with ex- 
to radial velocity variations (iMavor fc Quelozl[l995h . This treme caution, it is true that a large number of PCEBs 
finding generated large radial velocity surveys for exoplan- show a PP ar ent period variations that might be explained 
ets, and meanwhile more than 600 planets around main se- bv changes m the light travel time caused by a circumbi- 
quence (MS) stars have been identified using this method. narv § iant P lanet - We P resent an observational census of 

eclipsing, close-compact binaries and eclipse timing mea- 

Send offprint requests to: M. Zorotovic surements in the appendix. 
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If confirmed, PCEBs as hosts of giant planets are par- 
ticularly interesting because the binary host star passed 
through a special evolutionary phase: common -envelope 
(CE) evolution. As outlined by iPaczvnskl (ll976f) and de- 
scribed in much more detail by IWebbm k (1984) CE evo- 
lution occurs once the initially more massive star (the pri- 
mary) evolves. If it fills its Roche lobe during one of the gi- 
ant phases, when it has a deep convective envelope, dynam- 
ically unstable mass transfer to the less massive component 
(the secondary) begins. The secondary is not able to adjust 
its structure on the mass-transfer time scale. The mate- 
rial lost by the primary therefore initially fills the Roche 
lobe of the secondary and then grows to form a noncoro- 
tating CE that surrounds the core of the giant and the 
secondary star. Owing to drag forces within the envelope, 
orbital energy is extracted from the binary and transferred 
to the envelope, which dramatically reduces the separation 
between the core of the primary and the secondary star in 
a spiraling-in process, until the envelope becomes unbound 
and is ejected from the binary. The remaining system is 
a PCEB consisting of the core of the primary (a compact 
object) and an MS companion, in a close but detached or- 
bit. The short duration of the CE phase (< 10 3 yr) means 
that the mass of the secondary sta r is assumed to remain 
constant (|Hiellming fc Taanil Il9 91). and the mass of the 
compact object should be equal to the mass of the core of 
the giant at the onset of mass transfer (for more details of 
CE evolution see, e.g ., Ilben fc Liviol 119931: 1 Webbinkl [200l 
IZorotovic et~alll2010h . 

The special evolution of the host binary stars implies 
that the claimed circumbinary planets must have either 
survived the dramatic evolution of the host binary star 
or must have formed as a consequence of this evolution. 
Which of these two spectacular scenarios might have taken 
place in PCEBs with candidate circumbinary planets is cur- 
re ntly a completely open question (see, e.g., the discussion 
in lBeuermann et aLll2010D . Answering it based on hydrody- 
na mic simulations o f CE evolution, such as those presented 
bv lRicker fe Taaml (|2012l ). is currently impossible, because 
CE evolution is still relatively uncertain even without an 
embedded planet. 

In this paper, we present the results of binary popula- 
tion synthesis simulations and reconstruct the evolution- 
ary history of observed eclipsing PCEBs to characterize 
their MS+MS binary progenitors. Reviewing what is known 
about circumbinary planets and circumbinary planet for- 
mation around such MS+MS binaries, we conclude that 
the apparent period variations observed in PCEBs are very 
unlikely to be caused by giant planets that survived CE 
evolution, and we discuss alternative explanations. We start 
with a review of the observed sample of eclipsing PCEBs 
and the detected apparent period variations. 

2. Giant planets around compact binaries? 

In the past few years, it has been suggested that there are 
giant planets around several eclipsing, close-compact bina- 
ries to explain observed variations of eclipse timings. These 
discoveries have become possible thanks to both the in- 
crease in the number of known eclipsing systems in recent 
years and the improvement in the accuracy of the observed 

light curves and derived eclipse timings. 

Almost a decade ago, Schreiber & Giinsicke ()2003l ) an- 
alyzed 30 detached PCEBs when only 11 eclipsing systems, 



including the prototype V471 Tau, were known. The SDSS 
and intensive follow-up observations of WD+MS binaries, 
as well as the Catalina Realtime Transient Survey, signifi- 
cantly increased the number of known eclipsing, detached 
PCEBs. TableQ] lists the orbital parameters of 56 known 
eclipsing, detached PCEBs and separates systems with WD 
and sdB primary stars. We excluded eclipsing PCEBs that 
are in the center of a planetary nebula, because no accurate 
eclipse timings are available for those systems. A brief re- 
view of the observational history of the 56 eclipsing PCEBs 
can be found in the appendix. 

Currently, 13 eclipsing, detached PCEBs consisting of 
an sdB (or sdOB) primary and a MS or BD companion are 
known. For an amazing fraction of ~ 38 per cent, apparent 
period changes that might be caused by a third substel- 
lar body have been measured (systems in bold in Tablc[l]). 
Even more dramatic, as shown in the appendix, almost all 
systems that have been intensively followed up for more 
than ~ 5 years show period variations that might indicate 
the p resence of a thi rd body, with AA Dor the only excep- 
tion (|Kilkennvll20Tl . 

The sample of known eclipsing, detached PCEBs with a 
WD primary contains 43 systems. For four of them, third- 
body detections using eclipse timing measurements have 
been claimed. Very few eclipse timings have been published 
for almost all other systems, mostly because the systems 
have been discovered fairly recently. In addition, several 
early eclipse timings are not reliable. This means that, as 
in the case of sdB+MS eclipsing binaries, all systems with 
accurate eclipse timings (error < 10 sec) covering > 5 years 
show apparent period changes that may indicate there is a 
third circumbinary object. 

In the past few years, apparent period variations asso- 
ciated to the presence of a third body have also been ob- 
served for three cataclysmic variables (CVs): UZFor, HU 
Aqr, and DP Leo. CVs are PCEBs that evolved into a semi- 
detached configuration where a Roche lobe-filling MS star 
(or a BD) transfers mass on the WD. Eclip sing CVs wer e 
discovered more than a century ago (see, e.g.. lPogsonll 857). 
and up to now there are almost 200 eclipsing systems among 
the ~ 1000 CVs listed in the catalog of iRitter fc Kolbl 
(2003, V7. 16). A complete discussion of all these systems 
is beyond the scope of this paper, but for completeness, we 
list the parameters of the three CVs with claimed planet 
detections in Table[T] and briefly review the observational 
history of these three systems in the appendix. 

The orbital parameters obtained with the best fit avail- 
able for the claimed planets are listed in Table[5J It should 
be noticed that, at present, not a single set of orbital pa- 
rameters has been confirmed by new measurements. In ad- 
dition, some of th e claimed planetary systems seem to be 
unstable (see, e.g.. lHorner et al.ll201lHHinse et al.ll2012l) . 

Still, the large fraction of PCEBs that show apparent 
period variations is intriguing. Nine out of ten detached 
PCEBs, which corresponds to a fraction of 90 ± 9 per cent, 
with accurate eclipse timing measurements covering ~ 5 
years show clear variations that might indicate the presence 
of a third body. If these potential planets have formed prior 
to CE evolution, a similarly large fraction of the MS+MS 
progenitor binaries of PCEBs must also host giant planets. 
To evaluate whether this might indeed be the case, we need 
to derive the orbital and stellar parameters of the PCEB 
progenitors. 
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3. Binary population simulation 

To characterize the progenitors of PCEBs we performed 
binary population studies of PCEBs containing a WD or 
an sdB primary. We complement our simulations by recon- 
structing the evolutionary history of the observed sample 
using the code described in lZorotovic et al.l (|2011aT l. 

Binary population studies of WD+MS PCEBs have 
been previously perfo r med b y several author s (see, 
e.g.. Ide Kool fc Ritterl 119931: IWillems fc Kolbl 12004 



iDavis et al.l 12010ft . and 

have simulated the pop- 



iPolitano fc Weilerl I2006L 1200 
recently. IClausen et al.l (|201 
ulation of sdB+MS binaries. We here present the first 
population simulation of PCEBs including both systems 
with WD and sdB primary stars. In addition, it is the 
first time that the present-day population of PCEBs is 
simulated with a low value of the CE efficiency, including a 
fr action of the rec o mbina tion energy, based on the results 
of I Zorotovic et all (|2010l) . A full parameter study will be 
presented elsewhere. 

3.1. Initial conditions and assumptions 

We use a Monte Carlo code to generate an initial population 
of 10 8 MS+MS binaries. The mass of the primary star is 
di stributed acco r ding to the initial mass function (IMF) 
of iKroupa et al.l (|l993f ); i.e., the number of primaries with 
masses in the range dM\ is given by dN oc f(M\)dMi 
where /(Mi) is given by 



/(Mr 





0.29056M{ 
0.15571Mf 
0.15571Mf 



1.3 



Mx/Mq < 0.1, 
0.1 < Mi/M < 0.5, 
0.5 < Mi/M < 1.0, 
1.0 < Mi/M . 



(1) 



The mass of the secondary is assumed to be distributed 
according to a fiat initial-mass-ratio distribution, i.e. n(q) = 
constant, where q = M2/M1. The initial orbital separation 
a; follows the distribution 







^ (ai) " ^ 0.078636^ 



(2) 



ai/R Q < 3 or ai /R Q > 10 6 
3 < a-jR Q < 10 6 , 

(|Davis et al.ll2008f ). We also assign a "born time" (ibom) to 
all the systems, corresponding to the age the galaxy had 
when the system was born, and assume that each ibom 
between and the age of the galaxy (i ga i ~ 13.5 Gyr) is 
equally likely, which corresponds to assuming a constant 
star formation rate. 

Once the initial population has been generated, we use 
the l atest version o f the binary-star evolution (BSE) code 
from iHurlev et al.l (|2002f) to evolve the systems. We have 
slightly updated the BSE code: the critical mass ratio for 
dynamically stable mass transfer when the prima ry is in the 
Hertzsprung gap was changed from 4.0 to 3.2 (|Han et al.l 
l2003f ). and the rate of angular momentum loss due to mag- 
netic brakin g was multip l ied by the normalization factor 
provided by Davis et al.l (|2008r ). Finally, we corrected a 



small mistake found by Robert Izzard (private communi- 
cation)!]] 



1 Eq. 32 from IHurlev et alJ (|2002l 1. related to tidal effects in 
close binaries, was misspelled in the code. The correct equation 

is fconv — TTV%Tl 



fconv 



[>«]■ 



while in the code it was written 



The systems are evolved for t cvo \ = £ ga i — £bom, to obtain 
the current orbital and stellar parameters. We assume a CE 
efficie ncy of 0.25 based on the results of IZorotovic et al.l 
(|2010D and the same fraction of recombination energy is 
included to compute the binding energy of the envelope. 

3.2. Results 

After evolving the systems, we obtain more than 4 x 10 5 
detached PCEBs consisting of a WD and a MS star or 
BD secondary. The mass and core composition of the WD 
clearly depend on the evolutionary stage of the primary 
star when the CE phase occurs. 

If the primary fills its Roche lobe during the first giant 
branch (FGB), the core will be mainly composed of helium. 
After the envelope is expelled, the core will not be massive 
enough to ignite helium and will evolve into a low-mass 
(< 0.5M Q ), helium-core WD (He WD). 

If the CE phase occurs close to the tip of the FGB, 
the core of the primary may be massive enough to start 
helium burning following envelope ejection. The minimum 
core mass for which this o ccurs depends on the initial mass 
of the primary (see, e.g., lHan et all I2002D . When helium 
ignites after envelope ejection, the star becomes a hot sdB 
star, i.e. a helium-core-burning star with a very thin (< 
O.O2M ) hydrogen envelope and a typical mass of ~ 0.45 — 
O.5M . It remains in this stage for ~ 10 8 yr until helium is 
exhausted in the core. The low mass of the envelope means 
the star cannot ascend the asymptotic giant branch (AGB) 
and instead bec omes a WD (for a recent review of sdB stars 
see lHeberl[2TJ?51 . 

Finally if the primary fills the Roche lobe on the AGB 
the core of the primary will be more massive (> O.5M ), 
and it will become a carbon/oxygen-core WD (C/O WD) 
after the ejection of the envelope. Very massive primaries 
can also produce oxygen/neon-core WDs but for the sake 
of simplicity we here call all post- AGB WDs C/O WDs. 

The BSE code does not distinguish between He WDs 
and sdB stars; i.e., all the systems that enter the CE phase 
when the primary is on the FGB are assumed to evolve into 
PCEBs with He WD primaries. However, based on their 
current and initial masses, we can identify which of these 
systems are actually SdB stars: those with a mass higher 
than the minimum core mass (M™ m ) needed to ignite he- 
lium after the envelo pe is lost. To com pute M™ ln we here 
use the results from lHan et al.l ([20021 their Table 1), as- 
suming a Reimers' wind mass-loss law (with rj — 1/4) and 
solar metallicit}(3. 

We divide the PCEB systems into four types, depend- 
ing on the type of primary star they contain: 

i) He WD primary, i.e., systems where the primary filled 
the Roche lobe on the FGB and Af wd < M c min ; 

ii) sdB primary, i.e., systems where the primary filled the 
Roche lobe close to the tip of the FGB with M wd > M c min , 
which happened less than 10 8 yr ago; 

iii) post-sdB primary, i.e., like the latter but if the CE 
phase occurred more than 10 8 yr ago, which means that 
the primary has been a helium burning sdB but has al- 



2 lHan et all (|2002T ) only computed M c min for two fixed values 
of metaflicity (0.02 and 0.004). However, it can be seen from 
their Fig. 1 that the difference is not dramatic for stars with 
initial masses below ~ 1.8 M©, which is the case for almost all 
the systems we obtain (see our Fig[2]) 
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ready evolved into a WD; 

iv) C/O WD primary, i.e., systems that entered the CE 
phase when the primary was on the AGB. 
Owing to the short duration of the sdB phase, we expect 
to find very few systems with primaries in this stage. 

It should be noticed that sdB stars can also be formed 
in binaries through stable mass transfer. However, their 
orbital periods are expected to be typically about two 
orders of magnitude larg er than those of PCEBs (e.g., 
IPodsiadlowski et alJ[2008T) . so we do not take into account 
this formation scenario here. 



3.2.1. The predicted current PCEB population 

Table[3] summarizes the number of systems that are found 
in each of the four groups just described, as well as their av- 
erage orbital parameters. The subscript "f corresponds to 
the final (current) value. As expected, very few (~ 2%) 
systems with sdB primaries are still in this stage. The 
sample is dominated by systems with C/O WD primaries. 
The relative number of systems with each type of primary 
depends on the value assumed for the CE efficiency, but 
the shape of the distributions are not significantly affecte d 
(see, e.g.. iWillems fc Kolbll200l iPolitano fc Weilerf 12001 . 
The resulting orbital period, primary mass, and secondary 
mass distributions of the predicted current PCEB pop- 
ulation are shown in Fig.[T] The complete population of 
detached PCEBs is represented by the gray shaded his- 
tograms. Distributions separating different types of pri- 
maries stars are shown as color histograms. Since we have 
obtained a very low percentage of systems that still con- 
tain an sdB primary, we multiplied the number of systems 
in each bin by ten (in green). 

Systems in which the primary fills the Roche lobe at 
a more advanced evolutionary stage, which usually corre- 
sponds to a longer initial period, may also end up with 
a longer final period. This is clearly reflected in the av- 
erage final periods listed in Tabled where the larger or- 
bital periods are for systems containing C/O WDs, and is 
in perfect agreement wi th observations of WD+MS PCEBs 
(jZorotovic et al.ll2011bl ). However, it should be noticed that 
the the mean value (P/) of the period distribution for 
PCEBs containing C/O WDs is strongly affected by a huge 
tail towards very long orbital periods. This tail corresponds 
to systems that fill the Roche lobe very late on the AGB 
and whose orbital separation decreases very slightly during 
CE evolution owing to the large amount of recombination 
energy available within the envelope. This is also reflected 
in the high value of the standard deviation for the mean 
value. Therefore, we also give in Tabled the value of the 
median (P/), which is a better representation of the peak 
of the distribution and which still agrees with the results of 
(|Zorotovic et alj|2011bl) . 

As expected, systems containing C/O WDs have M\ t t > 
0.5M Q , while progenitors that filled their Roche lobe on the 
FGB have M\j < O.5M separated by a small gap. The 
lower boundary of this gap is given by the maximum core 
mass that a giant star can have at the end of the FGB (~ 
0.48M Q ). The upper boundary appears because the radius 
of the star at the beginning of the AGB is smaller than at 
the end of the FGB. The minimum WD mass for a star 
that filled its Roche lobe during the AGB is given by the 
core mass at which the radius first exceeds the maximum 



0.03 




-1 1 

log(P,[d]) 



0.6 0.9 
M„[M B 



Figure 1. Current orbital period {left), primary mass [mid- 
dle), and secondary mass {right) distribution of the simu- 
lated PCEB population, separated according to the type 
of primary star they contain. The gray shaded histograms 
represent the entire population. 



FGB radius. This happens for core masses ~ 0.51M(rJ3- The 
mass distributions of He WDs and sdB (or post-sdB) stars 
overlap slightly, because the minimum core mass to ignite 
helium depends on the initial mass, but He WDs are in 
general slightly less massive than sdBs. 

The distribution of companion masses is similar for all 
types of PCEBs, which is also consistent with the observa- 
tions. There are more and more systems with increasing sec- 
ondary masses, with a steep decline at M2 ~ 0.35Mq, i.e. at 
the boundary for fully convective secondaries. According to 
the disrupted magnetic braking theory, systems with more 
massive secondaries suffer from angular momentum loss due 
to magnetic braking, in addition to gravitational radiation 
and hence spend less time in the detac hed phase. This ef- 
fect h as been predicted theoretically bv lPolitano fc Weiler 



(2007) and observationally confirmed by ISchreiber et al 



As expected, the orbital parameters for systems 
with sdB and post-sdB primaries are similar, with systems 
with a post-sdB primary having slightly shorter orbital pe- 
riods since they have had more time to evolve towards 
shorter periods after the CE phase. 

3.2.2. Relating the initial separation and final system 
parameters 

In the context of detecting circumbinary planets around 
PCEBs, the initial parameters of the systems, and their 
relation with the current parameters is crucial. 

One of the basic assumptions for the CE phase is that 
the mass of the secondary star remains constant, there- 
fore the initial secondary masses are almost identical to the 
current secondary masses. The distributions of initial pri- 
mary masses and initial separations are shown in Fig. [31 
and the mean values are listed in TableUJ Again, systems 
with sdB and post-sdB primaries cover the same range of 
initial parameters. As in Fig JT] the number of systems with 
sdB primaries was multiplied by ten in Fig. [5] 



These limits slightly depend on the metallicity. 
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Table 3. Current orbital properties of the simulated PCEB population separated by the type of primary star. The last 
column gives the median value for Pf, while all the others correspond to the means (see text for details). 













'/ = («f) 






Primary 


N 


Percent 


Mij 


Mi 


Pf 


Pf 




xlO 4 


% 


[M a ] 


[M ] 


[d] 


[d] 


He WD 


11.79 


26.06 


0.41 ± 0.02 


0.65 ± 0.29 


1.57 ± 0.71 


1.03 ± 0.86 


0.78 


sdB 


0.14 


0.31 


0.46 ± 0.01 


0.65 ± 0.39 


1.42 ± 0.86 


3.23 ± 3.42 


2.01 


post-sdB 


6.87 


15.18 


0.46 ± 0.01 


0.61 ± 0.31 


1.34 ± 0.69 


1.81 ± 1.69 


1.32 


C/O WD 


26.45 


58.45 


0.59 ± 0.12 


0.67 ± 0.42 


1.15 ± 0.69 


50.75 ± 154.94 


2.54 



Table 4. Initial primary mass, mass ratio, and orbital separation of the simulated PCEBs with the different primary 
types previously defined. 



Primary 


Mi,, 
[Mq] 


* = te) 


a, 

[Rq] 


a, 
[AU] 


He WD 


1.20 ± 0.15 


0.53 ± 0.23 


238 ± 53 


1.11 ± 0.25 


sdB 


1.32 ± 0.23 


0.48 ± 0.26 


361 ± 51 


1.68 ± 0.24 


post-sdB 


1.39 ± 0.19 


0.45 ± 0.22 


358 ± 45 


1.67 ± 0.21 


C/O WD 


2.08 ± 0.93 


0.34 ± 0.21 


699 ± 356 


3.25 ± 1.65 




18 3 4 300 600 900 1200 

M,,,[MJ a,[R s ] 



Figure 2. Initial mass distribution of the primary (left) and 
initial separation (right) for the different types of PCEBs 
obtained with our simulation. The short vertical lines in the 
upper part of each panel correspond to the reconstructed 
values for the systems in Table[T] The dashed vertical line 
corresponds to 1 AU. 



All the primaries that fill their Roche lobe on the FGB, 
i.e. the progenitors of PCEBs with sdB and He-core WD 
primaries, descend from low-mass stars (M^; < 1.8M Q ), 
while the progenitors of C/O WDs are generally more mas- 
sive and cover a wider range of masses. 

Apparently, most PCEBs have formed from binaries 
with relatively small initial separations ai<1000i?©, i.e. 
< 5 AU. The dashed vertical line in the right-hand panel 
of Fig. [5] represents 1 AU for comparison. PCEBs contain- 
ing He WD primaries descend from very close binaries with 
at = 100 — 400i? Q , the progenitors of PCEBs with an sdB 
primary had initial separations of 200 — 500i?o, and PCEBs 
with C/O WD primaries formed mostly from systems with 



initial separations >5OOi?0. This is not surprising, because 
the radius of a giant star increases as the star evolves and 
the mass of the core grows, and therefore, larger separa- 
tions are needed for the Roche lobe to be filled with a more 
massive core. 

The relation between the initial binary separation and 
the initial and final stellar masses is shown in Fig.[3J There 
seems to be a relation with the initial primary mass and, 
more evidently, with the final WD mass. This is a conse- 
quence of the radius of a giant star depending almost ex- 
clusively on the core mass, i.e. on the mass of the future 
WD. This implies that the initial separation should also 
be larger in order to allow the star to evolve to this stage 
without filling the Roche lobe. 

To compare the predictions of our simulations with the 
observed sample of eclipsing PCE Bs, we use th e recon - 
struction algorithm presented in IZorotovic et al.l (|2011af) 
to obtain the initial orbital parameters for all the sys- 
tems in Table[T] with available periods and masses. The 
resulting values are included in Fig. [5] in the upper part 
of each paneQ. The sample of observed eclipsing PCEBs is 
naturally biased towards short orbital periods, hence also 
slightly towards small initial binary separations of just a 
few AU, but otherwise the properties of the simulated and 
reconstructed PCEBs are similar. 

In summary, we found that PCEBs in general and the 
observed population of eclipsing PCEBs in particular de- 
scend from binaries with primary masses of Mi ~ 1 — 3M Q , 
secondary masses from 0.1 — I.OMq and initial separations 
of at < 5 AU. 



4 We can only calculate the initial separations at the onset of 
the CE phase, which might differ from the actual initial separa- 
tion due to, e.g., mass loss and tidal effects that may affect the 
orbital separation prior to CE evolution. 
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Figure 3. Relation between the initial orbital separation and: initial primary mass (left), secondary mass (middle), and 
current mass of the primary (right). 



4. Planet formation and planets around MS+MS 
binaries 

In the previous sections we have derived the properties of 
the MS+MS binary progenitors of PCEBs. In what follows 
we review what is known about planet formation and cir- 
cumbinary planets around such MS+MS binaries. We start 
with a brief review of planet formation theories. 



4.1. Giant planet formation 

Giant planet formation theory has been developed early 
for the solar system and the pioneering works of the two 
still today comp e ting models, i.e. the core- accretion model 
(|Safronovl 119691: iGoldreich fc Ward! ll978Tl and the disk- 
instability model (|Cameronlll978f ). were published well be- 
fore the first exoplanets were identified. 

According to the disk-instability scenario, massive pro- 
toplanetary disks may fragment into dense cores that can 
contract further and form giant planets. G ravitational in- 
stabilities occur if Toomre's Q parameter (|Toomre| 1 19641 ) 
approaches unity and if the local cooling time is shorter 
than or equal to the orbital period at a given radius. 
Current models almost concor dantly (to the best of our 
knowledge, the only exception is Boss 2012) show that such 
conditions might be present in massive young (< 1 Myr) 
disks and on ly at distances from the c entral star > 20 — 
100 AU (e.g.. [B^Ievl[2009llCai et al.l[2010h . However, migra- 
tion of the planet and scattering can redistribute the plan- 
etary orbits, and one can therefore not exclude relatively 
close planets also being formed due to disk instabilities. 

According to the alternative scenario, the core-accretion 
model, larger protoplanetary cores are built by the accre- 
tion of planetesimals until the core is massive enough to 



accrete gas from the surrounding protoplanetary disk. The 
accretion of planetesimals is relatively slow, and the forma- 
tion of giant planets according to the core-accretion m odels 
takes around four million years (jAlibert et al.ll2005l) . The 
strengths of the core-accretion model ar e that it naturally 
expla ins the planet-metallicity relation (Fi scher fc Valentl 
2005), that gas giants formed by core accretion are enriched 
due to the accretion of planetesim als as was bee n found to 
be the case in our solar system (|Yound |2003|) , and that 
rocky planets can form due to basically the same mecha- 
nism. Furthermore, population models based on the core- 
accretion scenario predict mass and semi- major axis distri- 
butions in agreem ent with observations (|Ida fc Linl l2004t 
lAlibert et al] l2005l. Only very few mass ve planets around 
low-metallicity stars might be di fficult to be explaine d 
within the core-accretion scenario (Mordasi ni et al.l l2012'). 
The core-accretion model is currently the favored forma- 
tion model for most giant planets although giant planet 
formation due to disk instabilities cannot be excluded. 



4.2. Multiplicity and exoplanets 



Already almost two decades ago, iGhez et all (|l993tl were 
the first to find a difference between the binary separation 
distributions of classical accreting TTauri stars (CTTS) 
and disk-less weak- lined T Tauri stars (WTTS), in the sense 
that the fraction of WTTS to CTTS is increased for short 
projected separations < 50 AU. This finding has meanwhile 
been impressively confirmed by several studies using the un- 
precedented datab ase provided by the Spitzer survey (e.g., 
ICieza et al.l l2009[). and similar r esults are found for the 
transition disk (jcieza et al. 20121) and debris disk systems 
(jRodriguez fc Zuckermanll2012f) . The lack of binary stars 
with projected binary separations ~ 1 — 100AU among 
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samples of young stellar objects still hosting a disk may in- 
dicate a short lifetime of circumbinary disks around these 
close binaries. 

Indeed, very recently, iKraus et al.l (|201ll l2012h have 
presented results of a high-resolution imaging study of the 
Taurus Auriga star-forming region of mostly solar type 
stars (M ~ 0.25 — 2.5M ) and compared the fraction of 
pre-MS stars hosting disks at different ages. They clearly 
show that the paucity of binaries among CTTS, transition 
disks, and debris disks is explained by the short lifetime 
of most circumbinary protoplanetary disks, i.e. < 1 Myr. As 
this is short compared to the time scale of giant planet for- 
mation according to the core-accretion model, only a small 
fraction of these binaries should host circumbinary plan- 
ets if core accretion is indeed the main channel for forming 
giant planets. 

In addition to the short lifetimes of circumbinary disks 
around solar- type binary stars with projected separations 
> 1 AU, the search of planets around MS+MS binary stars 
within the field of NASA's Kepler mission revealed a low 
frequency of giant planets around even closer solar-type 
MS+MS binaries. The first pl anet around a pai r of MS 
stars, Keple r- 16, w as found bv iDovle et ail (|2011f ). Later, 
IWelsh et alj (|2012h found two more planet transits-like 
events around the MS+MS binaries Kepler-34 and Kepler- 
35, in a survey of 750 eclipsing binaries with periods rang- 
ing from 0.9 to 276 days, from Kepler's survey of MS+MS 
binaries. These three systems have binary periods longer 
than 20 days, while most of the systems observed by Kepler 
(> 80%) have shorter orbital periods. Considering that it 
is more likely to detect planets around short-period sys- 
tems, this suggests a much lower rate of gas giant planets 
around t he shorter per i od bi naries. Based on geometric ar- 
guments, IWelsh et al.l (|2012h infer a lower limit of - 1% 
of circumbinary giant planets with periods around 100-200 
days in close binaries (P = 0.9-50 days). 

In summary, the low frequency of circumbinary planets 
complements the apparently short lifetimes of circumbinary 
disks for projected separations of > 1AU. We conclude that 
evidence for small numbers of circumbinary giant planets 
around close solar-type binary stars is growing. Such binary 
stars are the progenitors of PCEBs, so that the high fre- 
quency of apparent period variations observed in eclipsing 
PCEBs (see Sect. [5]) is unlikely to be caused by giant plan- 
ets that survived CE evolution. In fact, the only remain- 
ing possibility to explain the period variations observed in 
nearly all PCEBs due to such first-generation planets re- 
quires that nearly all first-generation circumbinary disks 
form giant planets at large separations that are not yet de- 
tectable with Kepler. Furthermore, these planets must form 
due to the fast disk-instability mechanism that is not in con- 
flict with the short disk lifetimes mentioned above. While 
current observations cannot entirely exclude this possibil- 
ity, it appears unlikely given that a flat extrapolation of 
the observed planet separation distribution to separations 
of 3 — 10 AU around single stars predicts that only ~ 12 
per cent may host giant planets (jMarcv et al.ll2005f i. 

5. Discussion 

In the previous sections we obtained the following results 

— From the 56 eclipsing, detached PCEBs with MS or BD 
secondaries, accurate eclipse time measurements cover- 



ing more than around five years have been obtained for 
ten systems. All of these but A A Dor show apparent 
period changes that might indicate the presence of cir- 
cumbinary planets. 

— The progenitors of PCEBs are MS+MS binaries with 
initial se parations of ~ 0.5 — 5 AU, similar to tho se 
tested bv lKraus et~ai1 (|2012h and lWelsh et al.l (j2012ft . 

— The reduced lifetime of circumbinary protoplanetary 
disks does not allow the frequent formation of cir- 
cumbinary planets around close MS+MS binaries 
(jKraus et al.1 [2012) within the core-accretion model. 
This has been shown for projected separations ~ 1 — 
40 AU and is likel y to also be true for closer binaries 
(|Welsh et alJl20H . 

These results suggest that the apparent period changes 
in PCEBs are not likely to be caused by circumbinary gi- 
ant planets that formed in the circumbinary protoplane- 
tary disk around the initial MS+MS binary. Instead, if 
the variations are indeed caused by circumbinary plan- 
ets, the planets may have formed following CE evolu- 
tion in a second-generation disk. The idea of second- 
generation planets goes back to the first exoplanet dis- 
covered around the pulsar PSRB1257+12. According to 
iTavani fc Brookshawl (|l992t) , a circumbinary protoplane- 
tary disk may have formed around the pulsar from evapo- 
rated material that did not leave the system. More recent 
models start from supernova fallback disks and can repro- 
duce the observed eccentricities and masses of the plane - 
tary system around PSRB1257+12 (jHansen et al.l l2009h . 
It is therefore reasonable to suggest that also the claimed 
planets around PCEBs have a second-generation origin. In 
such a scenario, a protoplanetary disk must have formed 
following the CE ejection. This idea is supported by recent 
simulations of CE evolution indicating that envelope ejec- 
tion can be incomplete, i.e. up to ~ 10 per cent of the 
ejected material may remain bound to the PCEB and is 
likely to form a circ umbinary disk ()Kashi fc Sokerl 120111 : 
iRicker Sz Taaml [2012T ). If the apparent period changes in 
PCEBs are indeed caused by planets formed in such a 
second-generation protoplanetary disk, the planet forma- 
tion process must proceed on a short time scale since the 
cooling age of the most convincing PCEB with planets, NN 
Ser, is just a million years. It must also be very efficient be- 
cause nearly all PCEBs with accurate eclipse timings show 
apparent period variations. Below we discuss several op- 
tions of second-generation planet formation and add a note 
of caution concerning the third-body interpretation. 

5.1. Second-generation disk instabilities 

The disk-instability model offers planet formation on short 
time scales that could solve the disagreement with the 
circumbinary planets around very young PCEBs such as 
NN Ser. However, in general, the disk-instability model still 
has to ove r come serious problems, as recently discussed by 
IZhu et all (|2012h . The orbital separations of the claimed 
planet detections around PCEBs disagree with the predic- 
tions of the disk instability model. Almost all the claimed 
planets are in orbits with a<10 AU (see Tabled, while 
the conditions for forming planets according to the disk- 
instability scenario occur at distances > 20 AU from the 
central star. Significant planet migration and/or scattering 
would thus be necessary. Furthermore, given the high frac- 
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tion of PCEBs showing period variations would require that 
nearly all post-CE circumbinary disk are unstable and form 
giant planets. Although this scenario cannot be excluded, 
it is not clear why second-generation disks should be al- 
most always unstable, fragment, and form giant planets, 
while disks around young single stars and MS+MS binaries 
seem to form giant planets in only about ten per cent of 
all cases. Therefore, second-generation disk instabilities do 
not appear to be the most promising solution. 

5.2. Second-generation rapid core accretion 

Planet formation according to the core- accretion model is 
more likely for higher metallicities (e.g. iFischer fc Valentl 
2005), and therefore the production of massive elements 
during the AGB evolution might favor the formation of 
second-generation giant planets. Indeed, dusty disks sur- 
rounding binary post-AGB s tars are a well known phe- 
nomenon. IWaters et all ()l998h were the first to claim that 
such disks are similar to protoplanetary disks and may lead 
to the formation of planets. More recently, evidence has 
grown that dust-processing in these disks and in proto- 
planetary disks aro und young stars is probably very similar 
(jGielen et al.ll20(M) . Moreover, a debris second-generation 
disk around an is olated neutron star has been detected 
(jWang et al.l l20061 . Furthermore, the second-generation gi- 
ant planet formation process might be fast if the envelope 
that forms the second-generation disk is rich of heavy el- 
ements, especially C and O, th at trigger enhanced dus t 
production and planet formation (|Petigura fc Marcv 2 0111) . 
Significantly enhanced dust-to-gas mass ratios may be able 
to explain both the high frequency of circumbinary planets 
and the short time scale of the second-generation planet for- 
mation required to explain the youth of the WD in NN Ser 
(~ lMyr). Therefore, the second-generation, rapid core- 
accretion scenario seems to be a reasonable option to ex- 
plain the observed apparent period variations. 

Even better, based on the above, the second-generation 
disk model makes predictions that can be tested observa- 
tionally. Enrichment of the envelope with heavy elements 
is well established as occurring at the end of the AGB 
but should be completely absent on the early FGB. As a 
result, second-generation disks around PCEBs with high- 
mass C/O WD primaries should have much higher metallic- 
ities than around systems with low-mass He WD primaries 
and should therefore be found to more frequently host giant 
planet sQ. 

A borderline case is PCEBs with sdB primary stars. 
Some authors claim that dust pro duction may also occur 
very close to the tip of the FGB ()Bover et al.l I2010T) . pre- 
cisely where the progenitors of sdBs fill their Roche lobes. 
The evolution towards sdB stars is far from being com- 
pletely understood. Late He-flashes may well lead to mix- 
ing in systems that experienced significant mass loss, but 
details are missing. We know neither the details of CE evo- 
lution nor pulsations or late He-flashes that may occur at 
the tip of the FGB. All we know for certain is that sdB stars 



5 However, this prediction is only based on metallicity, which 
might not be the only difference in the formation of C/O and 
He-core WDs that could affect potential planet formation in a 
second-generation circumbinary disk. Other factors, such as the 
expelled mass that is higher for AGB progenitors, may compli- 
cate the issue and dilute the prediction. 



are burning He, which produces C/O, and one may specu- 
late that this leads to enhanced metallicities in a potential 
second-generation disk around PCEBs with sdB primary 
stars. 

The current data is inconclusive with respect to possible 
relations between the nature of the compact object and oc- 
currence of circumbinary planets. There are nine PCEBs 
and three CVs with period variations attributed to cir- 
cumbinary planets (see Table[T]). Among the PCEBs, five 
of them contain an sdB primary, three a C/O WD, and 
only one a WD with Af wc i < 0.5. For the last, RRCae, ap- 
parent period variations have been only recently detected, 
and still need to be confirmed. Also, the mass of the WD 
is high enough to be either a He WD or a post-sdB star. 
In the case of CVs, all of them contain C/O WDs, which is 
not sur prising given the observ ed mass distribution of WDs 
in CVs (jZorotovic et al.|[2011al) . These fractions agree with 
the second-generation disk scenario but are still far from 
providing clear constraints. Intense monitoring of a large 
sample of eclipsing PCEBs is required to test the predic- 
tions of second-generation planet formation. 

5.3. A hybrid first- and second-generation scenario 

Neptune-mass planet s seem to be more fr equent than giant 
Jupiter- like planets (|Borucki et al.|[20TTl ). at least around 
single stars. If this can be further confirmed, especially for 
planets around close MS+MS binary stars, it might be that 
second-generation planet formation is efficient because the 
remnants of lower mass planets survived CE evolution and 
serve as seeds for more massive planets. As these planets 
already have a few Earth masses, the formation of a giant 
planet might be faster than for first-generation giant plan- 
ets, and the short cooling age of the WD in some PCEBs is 
no longer a problem. The predictions of this scenario are a 
high fraction of low-mass planets around MS+MS binaries, 
which subsequently can grow due to the accretion of gas 
left by the CE phase. If this is the case, there should be 
no relation between the type of compact object in PCEBs 
and the frequency of circumbinary planets. While low-mass 
circumbinary planets a round MS+MS bin aries have indeed 
recently been detected (jOrosz et al.ll2012aHb1 ). proper statis- 
tics of the frequency of low-mass circumbinary planets are 
not yet available. 

5.4. A note of caution: are we really detecting planets? 

As outlined in the introduction, pulsar timings led to de- 
tecting the first confirmed exoplanet orbiting the pulsa r 
PSRB1257+12 dWolszczan fe Fraill[l99l [Wolszczanlll99l . 
so planets around compact binary stars do exist. However, 
one also has to kee p in mind that th e possible compan- 
ion to PSR 1829-10 tiBailes et al.lll991l) was retracted later 
dLvne fc Baileslll99i . 

The situation for PCEBs, and especially CVs, is similar. 
All we see are eclipse timing variations that can be repro- 
duced well by assuming the existence of an orbiting cir- 
cumbinary third body, but not a single set of eclipse times 
that confirms the parameters of a previously claimed plan- 
etary systems has been published so far. In contrast , some 
suggested planetary systems (e.g.. iQian et al.ll2010bL|201ltl 
turned out to be d ynamically unstable ( Horner et alj |2011; 
iHinse et al.l 120121 ) or drastically disagreed with more re- 
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cent high-precision eclipse timings (Parsons et alJ 1201 0b). 
It therefore remains an open question whether we indeed 
identified a new, large, and exciting population of circumbi- 
nary extrasolar planets or if perhaps an alternative process, 
such as the freq uently proposed Applegate's mechanism 
(|Applegatd I1992D or a so far unknown process acting in 
deeply convective secondary stars might be responsible for 
at least some of the observed timing variations. In this con- 
text, AA Dor, a PCEB that may contain a BD companion, 
is extremely interesting because it is so far the only PCEB 
with continuous high-precision eclipse time measurements 
that does not sh ow any signs of apparent period variations 
(|Kilkennvll201l|) . One might therefore speculate that indeed 
the apparent period changes are somehow related to the 
existence of the convective secondary present in all PCEBs 
and CVs with claimed planet detections instead of being 
caused by circumbinary planets. Observations of eclipsing 
PCEBs with other types of secondary stars, i.e. with a WD 
or a BD companion (like e.g., SDSSJ0820+0008), may give 
new insight into this alternative possibility. 

6. Conclusion 

By combining binary population models with recent ob- 
servational and theoretical results for the formation of cir- 
cumbinary giant planets, we have shown that the appar- 
ent period variations seen in virtually all close-compact 
binaries with good coverage (>5yr of mid-eclipse times 
with good accuracy) are unlikely to be explained by first- 
generation planets. The lifetimes of protoplanetary disks 
around MS+MS binaries are simply too sho rt to form gi- 
ant planets in most cases (jKraus et al.ll2012T) . and the ob- 
served frac tion of planet s arou nd MS+MS stars further con- 
firms this ([ Welsh et al .112013) . The remaining options for 
explaining the observed period changes are either second- 
generation planet formation or alternative explanations 
that do not involve the existence of third and fourth bod- 
ies. We proposed observational experiments to test both 
hypothesis. 

First, if not caused by circumbinary planets but the ac- 
tive secondary, the period variations should not be detected 
in close WD binaries with a second WD component or with 
a BD secondary. Second, if caused by second-generation 
planets, a clear relation between dust production in the 
envelope of the compact object progenitor and planet fre- 
quency is expected; i.e., planets should be more frequent 
around compact binaries with C/O WDs than around those 
with a He-core WD primary. Continuous monitoring of all 
the eclipsing PCEBs listed in the appendix and also of all 
new systems that will be discovered will shed light on the 
origin of the observed apparent period changes in PCEBs. 

Acknowledgements. MZ acknowledges support from 

comite mixto, Gemini/Conicyt (grant 32100026) and 
CONICYT/FONDECYT/POSTDOCTORADO/3130559. MRS 
thanks FONDECYT (project 1100782) and the Millennium Science 
Initiative, Chilean Ministry of Economy, Nucleus P10-022-F. We 
thank Steven G. Parsons for helpful discussions. 



10 



Zorotovic & Schreiber: Origin of apparent period variations in PCEBs 



Table 1. Orbital parameters for the currently known eclipsing PCEBs. Systems with suspected planets are in bold. We distinguish 
between detached systems containing sdB primaries, detached systems with WD primaries, and CVs. For the case of CVs we only 
list those with suspected planets (see Sect. [5] for details). 



System 


Alt. Name 


orb 


Mi 


M 2 


References 






[d] 


[M ] 


[M ] 






Detached sdB+MS/BD PCEBs 


HW Vir 


2M J1244-0840 


0.11671955 


0.485 ±0.013 


0.142+0.004 


1,2 




HS 0705-1-6700 


2M .T0710+6655 

Mil J. « 9 V 3 1 _l_ \ 3 | V J \ J • 9 • 9 


095646625 


0.483 


0.134 


3,4 


HS 2231+2441 


2MJ2234+2456 


0.110588 


0.47: 


0.075: 


5 




NSVS 14256825 


2M J2020+0437 


0.1103741 


-0.46 


-0.21 


6 




NY Vir 


PG 1336-018 


0.101015967 


0.459+0.005 


0.122+0.001 


7, S 


1,9 


2M1938+4603 


NSVS 05629361 


0.1257653 


0.48+0.03 


0.12+0.01 


10 




NSVS 07826247 


CSS06833 


0.16177042 


0.376+0.055 


0.113+0.017 


11 




BUL-SC16 335 


2MJ1809-2641 


0.12505028 


0.5: 


0.16: 


12 




PG 1621+4737 


2MJ1622+4730 


0.075 






13 




SDSSJ0820+0008 


GSC 0196.0617 


0.097 


-0.25 


0.045+0.03 


14, 15 








-0.47 


0.068+0.03 


15 




ASAS 10232 


2MJ1023-3736 


0.13927 


0.461+0.051 


0.157+0.017 


16 




AA Dor 


LB 3459 


0.261539736 


0.471+0.005 


o.o788ig:88g 


17,18 


EC 10246-2707 




0.118507993 


0.45+0.17 


0.12+0.05 


19 




Detached WD+MS PCEBs 


NN Ser 


2M J1552+1254 


0.13008014 


0.535+0.012 


0.111+0.004 


20, 


21 


V471 Tau 


2M J0350+1714 


0.52118343 


0.84+0.05 


0.93+0.07 


22, 


23 


QS Vir 


EC 13471-1258 


0.1507575 


0.78+0.040 


0.430+0.040 


24 




RRCae 


2M J0421-4839 


0.30370363 


0.440+0.022 


0.183+0.013 


25 




DECvn 


RXJ1326.9+4532 


0.364139315 


n ri +0.06 
u - oi -0.02 


0.41+0.06 


26, 


27 


GKVir 


SDSSJ1415+0117 


0.344330833 


0.564+0.014 


0.116+0.003 


28 




RXJ2130.6+4710 


2MJ2130+4710 


0.52103562 


0.554+0.017 


0.555+0.023 


29 




SDSSJ0110+1326 


WD 0107+131 


0.332687 


0.47+0.02 


0.255-0.380 


30 




SDSSJ0303+0054 




0.1344377 


0.878-0.946 


0.224-0.282 


30 




SDSSJ0857+0342 


CSS03170 


0.06509654 


0.51+0.05 


0.09+0.01 


31 




SDSSJ1210+3347 




0.12448976 


0.415+0.010 


0.158+0.006 


32 








U .OOOO t 11 




n 07<j_i_n nno 


97 




SDSSJ1435+3733 




0.125631 


0.48-0.53 


0.190-0.246 


30 




SDSSJ1548+4057 




0.1855177 


0.614-0.678 


0.146-0.201 


30 




CSS06653 


SDSSJ1329+1230 


0.08096625 


0.350+0.081 




33, 


34 


CSS07125 


SDSSJ1410-0202 


0.363497 


0.470+0.055 


0.380+0.012 


34, 


35 


CSS080408 


SDSSJ1423+2409 


0.3820040 


0.410+0.024 


0.255+0.040 


34, 


35 


CSS080502 


SDSSJ0908+0604 


0.14943807 


0.370+0.018 


0.319+0.061 


33, 


34 


CSS09704 


SDSSJ2208-0115 


0.1565057 


0.37 




35 




CSS09797 


SDSSJ1456+1611 


0.229120 


0.370+0.016 


0.196+0.043 


34, 


35 


CSS21357 


SDSSJ1348+1834 


0.2484 


0.590+0.017 


0.319+0.061 


34, 


35 


CSS21616 


SDSSJ1325+2338 


0.1949589 






33 




CSS25601 


SDSSJ1244+1017 


0.227856 


0.400+0.026 


0,319+0.061 


34, 


35 


CSS38094 


SDSSJ0939+3258 


0.3309896 


0.520+0.026 


0.319+0.061 


33, 


34 


CSS40190 


SDSSJ0838+1914 


0.13011232 


0.390+0.035 


0.255+0.040 


33, 


34 


CSS41631 


SDSSJ0957+2342 


0.15087074 


0.430+0.025 


0.431+0.108 


33, 


34 


WD 1333+005 


SDSSJ1336+0017 


0.1219587 






33 




PTFEB11.441 


PTF1 J004546.0+415030.0 


0.3587 


0.51+0.09 


0.35+0.05 


36 




PTFEB28.235 


PTF1 J015256.6+384413.4 


0.3861 


0.65+0.11 


0.35+0.05 


36 




PTFEB28.852 


PTF1 J015524.7+373153.8 


0.4615 


0.52+0.05 


0.35+0.05 


36 




KIC-10544976 


USNO-B1.0 1377-0415424 


0.35046872 


0.61+0.04 


0.39+0.03 


37 




SDSS J0821+4559 




0.50909 


0.66+0.05 


0.431+0.108 


34, 


38 


SDSS J0927+3329 




2.30822 


0.59+0.05 


0.380+0.012 


34, 


38 


SDSS J0946+2030 




0.252861219 


0.62+0.10 


0.255+0.040 


34, 


38 


SDSS J0957+3001 




1.92612 


0.42+0.05 


0.380+0.012 


34, 


38 


SDSS J1021+1744 




0.14035907 


0.50+0.05 


0.319+0.061 


34, 


38 


SDSS J1028+0931 




0.23502576 


0.42+0.04 


0.380+0.012 


34, 


38 


SDSS J1057+1307 




0.1251621 


0.34+0.07 


0.255+0.040 


34, 


38 


SDSS J1223-0056 




0.09007 


0.45+0.06 


0.196+0.043 


34, 


38 


SDSS J1307+2156 




0.216322132 




0.319+0.061 


34, 


38 
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Table 1. Orbital parameters for the currently known eclipsing PCEBs. Systems with suspected planets are in bold. We distinguish 
between detached systems containing sdB primaries, detached systems with WD primaries, and CVs. For the case of CVs we only 
list those with suspected planets (see Sect. [5] for details). 



System 


Alt. Name 


orb 


Mi 


M 2 


References 






[d] 


[M ] 


[M ] 




SDSS J1408+2950 




0.1917902 


0.49±0.04 


0.255±0.040 


34, 38 


SDSS J1411+1028 




0.167509 


0.36±0.04 


0.380±0.012 


34, 38 


SDSS J2235+1428 




0.14445648 


0.45±0.06 


0.319±0.061 


34, 38 


CVs 


UZFor 


2MJ0335-2544 


0.08786542 


-0.71 


-0.14 


39, 40 


HU Aqr 


2M J2107-0517 


0.08682041 


0.80±0.04 


0.18±0.06 


41, 42 


DP Leo 


RX J2107.9-0518 


0.06236286 


1.2: 


0.14: 


43, 44 








0.6: 


0.09: 


45 



References. (1) Lee et all ([2009), (2) 


Bcuermann et al.l 


(2012a|). (5) 0stensen et al. (1200/1 (6) 


Wils et al. 


(2007 



(201 2M). (3 s ) iDrechsel et all (l200l. (4) 



9) iQian et all (l2012bD. ( 1 01 10stensen et al.l (120101)7 (11 



(20K 



Klepp fc Rauch| (l20ll . (19 



O'B rien et al.l (l2001h. ( 2 3) iKundra fc Hric 



(34) 
(38) 



(7) IVuckovic et al.l (l2007h. (8) 
jFor et al.l (l2010D. (12) i Polubek et a 



Beuermann et al.l 
Charpinet et all (120081)^ 



(14) iGeieretaLl (12012ft . (15) iGeier eTaLl (l2011b l). (16) ISchaffenroth et al.l (12011 



(18) 

(22) , _ 

(26) van den Besselaar et alj (120071) . 
(1200411 . (30) IPvrzas et al.l (12009ft 



Barlow et al 



iParsons et all (l2010al) 



27) iParsons et al 



lebassa-Manscrgas et al 
Parsons et al.l (l2012aD .~ (39) 



. (3 1) IParsons et al.l (1 20 
20121) . (35) iDrake et al ~ 



(20) 

( 24) lO'Donoghue et al.l (1200 
' 2010bD . (28) IParsons et al 



(32) IPvrzas et al 



(42) ISchwope et all (|201lH . (43) iPandel et al 
Notes: Very uncertain values are followed by 



Bailey fc Cropper 



201 Oh . (3 6) lLaw et al 
(11991. (40) iPotter et al.l 




1 (120071 ) (13) IGeier et all 



). (17) Kilkenny 



Bcuermann ct al 
(25) iMaxted et al 



( 2002), (44) iBeuermann et al.l (|2oTll ). (45) ISchwope et al.l (|2002). 




(I2012bh. (29) IMaxted etall 



Backhaus ct al 



(37) \M menara et al 



(41) ISchwarz et al 




It generally means that the mass was assumed and not derived. 



12 Zorotovic & Schreiber: Origin of apparent period variations in PCEBs 

Table 2. Best fits of the orbital parameters for the currently claimed planets around eclipsing PCEBs. 



Name 


Msin(i) 
[Mj] 


P 

H 


asin(i) 
[AU] 




e 




Ref. 


Notes 


HW Vir c 


14.3±1.0 


12.7±0.2 


4.69±0.06 





40±0 


10 


1 




HW Vir d 


30-120 


55±15 


12.8±0.2 




0.05: 




1 


* 


HS0705+6700 c 


31.5±1.0 


8.41±0.05 


3.52 





38±0 


05 


2 


* 


HS2231+2441 c 


13.94±2.20 


15.7 


~ 5.16 








3 




NSVS14256825 c 


2.8±0.3 


3.49±0.21 


1.9±0.3 





00±0 


08 


4 




NSVS14256825 d 


8.0±0.8 


6.86±0.25 


2.9±0.6 





52±0 


06 


4 




NY Vir c 


2.3±0.3 


7.9 


3.3±0.8 








5 




NY Vir d 


2.5: 


>15 


>5.08 








5 




NN Ser c 


6.91±0.54 


15.50±0.45 


5.38±0.20 




0.0 




6 




NN Ser d 


2.28±0.38 


7.75±0.35 


3.39±0.10 


o 


20±0 


02 


6 




V471 Tan c 

V J- I ±. lull 


46-111 


33 2±0 2 


~ 12 6-128 


o 


26±0 


02 


7 


* 


QS Vir c 


9.01 


14.4 


- 6.32 




0.62 




8 




QS Vir d 


56.59 


16.99 


- 7.15 




0.92 




8 


* 


RR Cae c 


4.2±0.4 


11.9±0.1 


5.3±0.6 









9 




UZ For c 


6.3±1.5 


16+3 


5.9±1.4 





04±0 


05 


10 




UZ For d 


7.7±1.2 


5.25±0.25 


2.8±0.5 





05±0 


05 


10 




HU Aqr c 


7.1 


9.00±0.05 


4.30 





13±0 


04 


11 




DP Leo c 


6.05±0.47 


28.01±2.00 


8.19±0.39 





39±0 


13 


12 





Refere nces. (I) IBeuermann et all (l2012bl). (2) iBeuermann et al.l (l2012al). (3\lQian et all (l2010cD. (A) lAlmeida et all 

f8l lAlmeida fc Jablonskil (l2011h 



(I2012D. | , 5)lQian et al.l (l2 012bl). (6) IBeuermann et al.l (l2010D. ^llKundra fc Hrid (l2011h. (8) lAlmeida , k Jablons 
(9) iQian et al.l (|2012al ) (lOl lPotter et al.l (|201lh . (11 ) iGozdziewski et al.l (|2012h . (12) IBeuermann et al.l (|201lH 



* The claimed third body is more consistent with a BD than with a planet. 
Notes: Very uncertain values are followed by ":" . 
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Appendix A: Notes on individual systems 

A.l. Eclipsing sdB+MS/BD PCEBs 

HW Vir is the prototype of HW Vir-like systems consist- 
ing of an sdB star with a late-type MS or BD companion. 
HWVir systems have short periods (~ 2 — 3 hr), and the 
light curves show very sharp primary and secondary minima 
an d a strong reflection effect . HW Vir itself was discovered 
bvlMenzies fc Marand (|1986( ), who determined the period. 
Kilkcnnv ct al. (1994) observed a change in the orbital pe- 
riod for the first time, w h ich motivated fu r ther s t udies (e.g., 
iCakirli fc DevIeD Il999l: iWood fc Safferl fl999l: iKiss et alJ 
1200(1: iKilkennv et aD l2000aL 120031: llbanoglu et aTll2004h . 
Many possible explanations have been proposed and dis- 
cussed, converging towards the existe nce of a th ir d ob - 
ject with a long period and low mass. iQian et aD (|2008l ) 
present new mid-eclipse times obtained from 2006 to 2008 
th at show some deviatio n fro m the sinusoida l fit pr oposed 
bv IKilkennv et all (|2003l ) and llbanoglu et all (|2004l ). They 
suggest a combination of cyclic variations plus a contin- 
uous decrease in the orbital per iod, which may reveal a 
fourth object with a long period. iLee et all (2009) present 
41 new mid-eclipse times taken from 2000 to 2008 and com- 
bine these with data from the literature. The Observed mi- 
nus Calculated (O-C) diagram of the orbital period span- 
ning more than 24 years shows a combination of two si- 
nusoidal variations, probably produced by the presence of 
two substcllar companions, plus a continuous period de- 
crease that is to o strong to be caus e d by gr avitational radi- 
ation. Recently, iBeuermann et all (|2012bl) have published 
26 new mid-eclipse times obtained between 2008 February 
and 2012 Febru ary, which deviate significantly from the 
I Lee et all (|2009h pr ediction. They a lso find that the so- 
lution presented bv ILee et al.l (|2009l ) is unstable and pro- 
pose a new solution involving two companions to HWVir: 
a planet and a BD or low-mass star. 

HS 0705+ 6700 is an sdB+dM binary di scovered to be 
a deta c hed short-period ecl ipsing system bv lDrechsel et al.l 
(|200lh . IQian et all ([2009b) obtained 38 mid-eclipse times 
between 2006 and 2008, updated the ephemeris, and 
performed an O-C diagram including 31 mid-eclipse 
times obtained from the literature since 2000 October 
(iDrechsel et alJ 120011; iNiarchos et~aD 120031: iNemeth et aD 
I2005t iKruspe et al.ll2007l ). Thev detected cyclic variations 
that were attributed to the light-trav el time effect pro - 
duced by the presence of a third object. [Qian et all (|2010d ) 
present new mid-eclipse times taken in 2009, propose the 
existence of a continuous decrease in the orbital period, 
and derive a mass for the third b ody corresponding to a 
substellar object, probably a BD. ICamurdan et alJ (|2012l ) 
obtained new mid-eclipse times in 2010 December, which 
disagree with a long-term period decrease. They used a si- 
nusoidal fit to adjust a third object in the system, which 
may be substellar or a very low-mass sta r , depen ding on the 
inclination. Recently, IBeuermann et all (|2012al) have pub- 
lished new mid-eclipse times obtained between 2009 August 
and 2011 December, updated the ephemeris, and derive the 
parameters for a possible third body that seem to be more 
consistent with a substellar object. 

HS 2231+2441 was discovered to be an eclipsing sdB 
with a low-mass, p robably substellar, companion by 
IQstensen et alJ (|2007f ). who determined the ephemeris and 
estimated the masses. Mid-eclipse times monitored since 



2006 (jQian et al.l l2010cJ) reveal a continuous decrease and 
cyclic variations in the orbital period. The cyclic oscilla- 
tion suggests ther e is a third object in the system, and 
IQian et all (|2010d ) estimate it may be a BD. 

NSVS 14 256825 wa s found to be an eclipsing sdB+dM 
binary by IWils et aD (|2007). who presented 19 pri- 
mary mid-eclipse times obtained between 2007 June and 
September and derived the ephemeris and orbital param- 
ete rs. The sy s tem ha s also been monitored since 2006 
bv IQian et alJ (|2010cD . whose preliminary results suggest 
a cyclic change in the orbital period. New mid-eclipse 
times obtained bet ween 2010 September and 2011 October 
were published by IKilkennv fc Koenl (|2012f) . They found 
an increase in period, but a long baseline is needed to 
see whether this behavio r is cyclic. In a parallel work, 
IBeuermann et all (|2012af ) present new mid-eclipse times 
obtained between 2009 July and 2011 October, which reveal 
an abrupt and continuous increase in the period after 2009. 
This is interpreted by the authors as the possible response 
to a t hird body in a highly elliptic orbit. In a very recent 
work. [Almeida et al.l (|2012 ) have reanalyzed the system in- 
cluding ten new mid-eclipse times between 2010 July and 
2012 August and find that the variations observed in the 
O-C diagram can be explained by two circumbinary giant 
planets. 

NY Vir is an s dB+dM binary whos e eclipsing na- 
ture was revealed by Kilkenny et aD (|1998l ). Subsequently, 
IKilkennv et alJ (|2000a ) present mid-eclipse times from 1996 
to 1999, determined the ephemeris, and fi nd no significant 
changes in the orbital period. However, Kilkenny] (|2011f) 
combine the previous results with new mid-eclipse times 
taken between 2001 and 2010 and find a continuous period 
decrease, which is too high to be due to gr avitational radia- 
tion. T he strong decrease was also seen bv ICamurdan et aD 
(|2012f ). who have presented new mid-eclipse times from 
2009 and 2011 and adj usted a downward parabola to the 
O-C diagram. Recently. IQian et aD (|2012bh have combined 
the previous results with nine new mid-eclipse times ob- 
tained in 2011 May, updated the ephemeris, and propose 
that the O-C diagram can be adjusted by a downward 
parabola plus a periodic variation produced by a planet 
orbiting the primary. They also suggest that the continu- 
ous period decrease may be part of a cyclic variation that 
may indicate the presence of a fourth object. 

2M 1 938+4603 was f o und t o be an eclipsing sdB+dM bi- 
nary bv lOstensen et all (|2010l ). who derived the ephemeris 
and orbital parameters. The available mid-eclipse times 
only cover the period between 2008 June and 2010 May 
and therefore no conclusion about apparent period varia- 
tions can be drawn yet. 

NSVS 07826247 is the longest period s dB+dM binary 
known to be eclipsing so far, discovered bv lKellev fc S haw 
(20071). Mid-eclipse times have been p ublished bvlFor et al. 
J2OI0I 2008 February to 2009 March ). iLiving fc Shengbang 
(|2010L 2009 March to August), and lBackhaus et alJ (|2oT2T 
2011 February to October). No evidence of period variation 
has been found so far, but more observations are needed in 
order to discard variations. 

BUL-SC16 335 was found to be an eclipsing system 
by iPolubek et al.1 (|2007f ). who suggest that it might be 
an HW Vir-like system based on the appearance of the 
light curve. They derived th e ephemeris and es t imate some 
of the orbital parameters. iTello fc Jablonsk i (2010) red- 
erived the orbital parameters and find disagreement with 
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iPolubek et al.l (|2007l ). The parameters are still poorly con- 
strained and no mid-eclipse times are available so far. 

PC 1621+4737 has been found in the course of the 
MUCHFUSS (massive unseen companions to hot faint un- 
de rluminous stars fro m SPSS. IGeier et all 1201 lah project 
bv lGeier et al.l (|2010f ). who determined the period and es- 
timated that the companion to the sdB star is a very late- 
type MS star or a BD. The system parameters are not con- 
strained and there are still no publishe d mid-eclipse times. 

SDSSJ0820+0008 was discovered bv lSchaffenroth et al.1 

(| 2 1 If) , who detected the short period radial velocity varia- 
tions, observed the eclipses in the light curve indicating an 
HW Vir-like systems, and determined the period. The mass 
of the sdB primary covers a wide range of possible solutions, 
and the co mpanion seems to b e a BD. The ephemeris was 
derived by IGeier et all (|2011bl ). Base d on high-re so lution 
spectroscopy with ESO-VLT/UVES IGeier et al.l (|2012h 
confirm the substellar nature of the companion. They also 
find a significant shift in the system velocity with respect to 
the previous study, which may be produced by a third ob- 
ject in the system. However, more obser vations are needed. 

AS AS 10232 was discovered by ISchaffenroth et al.l 
(|2011l) . who determined the period and observed the pri- 
mary eclipse in the light curve. The reflection effect, typi- 
cal of HW Vir-like systems, is clearly observed. The orbital 
parameters were derived based on photometry and spec- 
troscopy, and no mid-eclipse time has been published so 
far. 

A A Dor is a short- period eclipsing bina ry contain- 
ing an sdOB primary ([Kilkenny et all I1978D . Although 
many subsequent inv e stigat i ons have been publish e d (e.g., 
Kilkenny etall I1979L 119811: iKudritzki et al.l Il982t iRauch 
2000 1 : Iffilditch et alj|2q03t iFleig et al.ll20q8t IVuckovic et al. 
2008; iRucinskil 120091: iMuller et all 120101: Klepp fc Rauch 



20111) . the nature of the companion is still not clear, with 



some authors favoring a BD companion while others fa- 
vor a very low-mass M star. Eclipse timings were stud- 
ied by Kilkenny et al.l ()2000b|) , who found no e vidence of 
any p eriod variations between 1977 and 1999. Ki lkenny! 
(| 2 llh increased the baseline by including 13 new primary 
eclipse timings obtained between 2000 and 2010, updated 
the ephemeris, and confirmed the stability of the orbital 
period. 

EC 10246-2707 has recently been report ed as an eclips- 
ing sdB+dM binary by iBarlow et al.l (|2012f) . who estimate 
the orbital parameters, derive the ephemeris, and present 49 
mid-cclipse times covering 15 years between 1997 February 
and 2012 June. The O-C diagram reveals no secular changes 
larger than 10 -12 s s _1 in the period. However, the rela- 
tively low precision of the timings does not allow to rule 
out small-amplitude variations such as e.g. those observed 
in NN Ser. Additional observations of eclipses with high pre- 
cision are needed. 



A.2. Eclipsing WD+MS PCEBs 

NNSe r was discovered to have deep eclipses by lHaefnerl 
(|!989lL who classified the system as a DA+dM bi- 
nary and measured the orbital period. iBrinkworth et all 
(2006) performed high-time- resolution photometry with 
ULTRACAM and detected a decrease in the orbital pe- 
riod that Applegate's mechanism fails to explain. They 
suggest it may be related to the presence of a third body. 
Qian et al. (2009a) propose a sinusoidal fit to the O-C di- 



agram. Later, iParsons et al.l (|2010al lbT) presented the result 
of eclipse observations performed with ULTRACAM since 
2002, which disag ree with the sinusoidal fit proposed by 
Qian et all (|2009al ). However, they still consider that a third 
body may be the caus e of the observed period changes. 
iBeuermann et al.l (|2010f ) monitored the system during 2010 
and combined their new mid-eclipse times with all the pre- 
viously published times, reanalyzing those where it was nec- 
essary. They conclude that the large amplitude variations 
observed in the period can only be caused by a third body, 
and even suggest that the best model is obtained with two 
giant planets around the binary. However, the existence of 
the fourth body is still rather uncertain. 

V4 71 Tau is a DA+dK2 sy stem discovered to be eclips- 
ing bv lNelson fc Young] () 1970l L^ Period variations have been 



observed for a long time (e.g., Lohsenlfl974fl . and many pos- 
sibilities were suggeste d in the past to explain the shape of 
the O-C diagram (e.g., Ibanoglu et al. 1994 [ Ibanqglu et al.l 



120051 : lGuinan~fc Ribas 2001; Kamin ski et al'l l2007f. such as 
perturbations by a third body in a long-period orbit, apsi- 
dal motion due to a low orbital ec c entric ity, or even mass 
transfer. Recently, iKundra &; Hrid (|201lh have detected a 
change in the O-C diagram trend. This allows them to ex- 
clude mass transfer and other models that predict a further 
increase in the O-C value. After modeling the system, the 
authors find that the third component could be a BD or a 
very low-mass star, at a period of 33.2 years. 

QS Vir was discovered in the Edinburgh- Cape faint blue 
objec t survey of high galactic latitudes (Kilk enny et al.l 
1997 ), where the eclip ses revealed its binary nature. 



lO'Donoehue et all (|2003l ) sugge st that it is a hibern ating 
CV, which was questioned by Ribeiro et al.l |2010) and 
ruled out bv lParsons et al. who confirm that it is a 

detached system and not a hibernating CV based on high- 
resolution U VES spectr a . Orbi tal period variations were 
analyzed by lOian et all ()2010bl ) , who combined new and 
previously published mid-eclipse times, and propose that 
there is a giant planet a nd a continuous decrea se in period 
due to magnetic braking. IParsons et al.l ()2010b|) update the 
O-C diagram by including ULTRAC AM photom e try an d 
find strong disagreement with the fit of lQian et al.l ()2010bl ). 
They conclude that the decrease in orbital period is part of 
a cyclic variation that cannot be explained by Applegate's 
mechanism. A thir d body seems to be t h e mos t probable 
solution. Recently, lAlmeida fc Jablonskil (|2011|) have pre- 
sented new mid-eclipse times and suggest that the best fit 
in the O-C diagram is obtained with a model with two cir- 
cumbinary bodies, most likely a giant planet and a BD. 
RR Cae is a WD+dM binary discovered as a high proper 



1955). The eclipses were first an- 



motion obje ct by Luvtenl 

nounced by iKrzeminski (|1984|) . an d further obser v ations 
of t he eclipses we r e pre sented by iBruch fc Diazl (I199E 



and iMaxted et all (|2007l ). The later paper updated the 
ephemeris and finds no evidence of variati ons in the or- 
bital pe riod on a long time scale (~ 10 yr). IParsons et al.l 
(|2010bl) performed ULTRACAM photometry for the sys- 
tem, obtained two new mid-eclipse times, and combined 
these with all the previous eclipse times available in or- 
der to study possible period variation. They obtained a 
roughly sinusoidal variation in the O-C diag ram, which 
can be explained via Applegate's mechanism. iQian et all 
(2012aJ) have recently obtained six new mid-eclipse times 
that combined with those from the literature, show some 
evidence of a third object, a giant planet, and even possi- 
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ble evidence of a fourth companion. The last needs to be 
confirmed. 

DE CVn is a bright eclipsing WD+dM binary dis- 
cover ed as an X-ray source by R OSAT ([Voges et al.l 
|1999|) . Recently, [Parsons et all (|2010bf ) have obtained high- 
time-resolution photometry with ULTRACAM to obtain 
an accurate ephemeris. They combined their new mid- 
eclip se times with older times available in the litera- 
ture (iRobb fc Greimelll997t Ivan den Besselaar et al.1l2007t 
iTas et al.1 l200l " in order to study possible period varia- 
tions. However, only the ULTRACAM data are reliable, 
and there are still too few to analyze possible long-term 
period changes. 

GKVir was discovered by iGreen et al.l |lT)7->;. who 
listed nine mid-eclipse times from 1975 April to 1978 
February, and determined the orbital period. For more 
than twenty years there were no new eclipses observed. 
Between 2002 and 2007 IParsons et all (|2010bft observed 
seven primary eclipses with ULTRACAM and improved the 
ephemeris. After c ombining their points with those from 
IGreen et al.l (|1978l ), they observe d a period i n creas e and 
a slight variation in O-C times. iDrake et al.l ([2010) also 
published one mid-ec lipse time for 2005 April. Recently, 
IParsons et al.1 (|2012bl ) have obtained a new high-precision 
mid-eclipse time on April 2010 using ULTRACAM, which 
shows a clear deviation from linearity. The magnitude of 
the period change is small and therefore can be caused 
by Applegate's mechanism or due to a third body in the 
system. More data is needed before the true cause of this 
period change can be established. 

R X J 21 30. 6 +4 71 was discovered to have eclipses by 
iMaxted et al.l (|2004l ), who determined the ephemeris and 
published three ULTRACAM mid-eclipse times obtained 
in 2002 and 2003. No new mid-eclipse times are available. 

SDSSJ0110+1326 and SD SSJ0303+0054 were iden- 
tified as eclipsing binaries by iPvrzas et ail ([2009), who 
studied the eclipses between 2006 September and 2007 
October, derived the ephemeris, and listed four and seven 
mid-eclipse times, respectively. Subsequently, Parsons et al. 
(|2010bl) analyzed accurate eclipses obtained during 2007 
October using ULTRACAM and listed one and three new 
mid-eclipse times, respectively. T hey found some de via- 
tions from the ephemeris derived bv lPvrzas et al.l (|2009l ) for 
SDSSJ0110+1326, while the eclipses for SDSSJ0303+0054 
appear to be consistent in the two studies. However, 
given the short baseline and the large unc e rtaint y in the 
mid-eclipse times derived by IPvrzas et al.1 (|2009h further 
accurate observations are needed before long-term pe- 
riod c hanges can be explored. Recently, iBackhaus et al.1 
(|2012h have presented six new mid-eclipse times for 
SDSSJ0303+0054 between 2011 August and November, 
which increases the baseline for this system. The new 
eclipses are still consistent with linear ephemeris, however 
the authors still do not exclude long-term period variations. 

SDSSJ0857+0342 was first listed as an eclipsing system 
bv lDrake et al.1 (|2010D , who observed the regular eclipses as 
part of the Catalina Realtime Transient Survey and deter- 
mined the short orbital period, which makes it the closest 
detached WD+dM binary. T hey published a mid-e clipse 
time for 2005 April. Recently, IParsons et al.l (|2011al ) have 
presented nine new mid-eclipse times obtained between 
2010 December and 2011 January with ULTRACAM, and 
updated the ephemeri s that is so f ar con sistent with the 
eclipse time listed in IDrake et al.1 (|2010l) . Independently, 



IBackhaus et al.l ((2012) obtained seven more mid-eclipse 
times between 2010 November and 2011 October. Due to 
the uncertainty in Drake's eclipse, possible period changes 
cannot be discarded so far. 

SPSS 'J 12 1 0+33 J 7 was discovered to be eclipsing by 
IPvrzas et al.1 (|2012l ). who obtained high-time- resolution 
photometry for nine eclipses between 2009 April and 2011 
June using RISE on the Liverpool Telescope. They deter- 
mined the orbital period and the ephemeris, however the 
mid-eclipse times still cover a short period of observation 
to study changes in the period. 

SDSS J1212-0123 was listed as an eclip sing WD+dM 
binary bv lNebot Gomez-Moran et al.l ([20091 ) , who obtained 
six mid-eclipse times between 2007 January a nd 200 8 May 
and derived the ephemeris. IParsons" et al.1 (|2012bh used 
ULTRACAM to obtained the first high-precision mid- 
eclipse in 2010 A pril, which was consistent with th e mid- 
cclipse times from lNebot Gomez-Moran et al.l (|2009t ) . They 
updated the ephemeris using the new accurate eclipse. 

SD SSJ1435+3733 is a par tially eclipsing binary discov- 
ered bv lSteinfadt et ail (|2008f ), who observed three eclipses 
i n 200 7 May and June. In a parallel study IPvrzas et al.l 
(2009) independently identified the system, observed seven 
eclipses between 2007 February and May, and updated the 
ephemeris. 

SDSSJ1548+40 57 was found to be eclipsing by 
IPvrzas et al.1 (|2009f) . who derived the ephemeris and listed 
seven mid-eclipse tim es from May to 2008 July. Recently, 
Back haus et al.l (|2012tl presented six new mid-eclipse times 
between 2011 May and August, which are consistent with 
linear ephemeris but are still not enough to exclude long- 
term period variations. 

CSS06653, CSS080502, CSS21616, CSS38094, 
CSS40190, CSS41631, and WD 1333 +005 are all WD+ dM 
binaries discovered to be eclipsing bv lDrake et al.l (|2010l ) as 
part of the Catalina Realtime Transient Survey. They only 
list one mid-ecl ip se tim e in 2005 April or May. Recently, 
IBackhaus et~aT1 (|2012h present new mid-eclipse times 
obtained between 2011 January and October for all these 
systems, and derived the ephemeris again. The accuracy of 
Drake's eclipses is not good enough to confirm or exclude 
long-term period variations. However, most of the systems 
seem to be consistent with linear ephemeris (within a 2-er 
error), except for CSS06653 and WD 1333+005, which 
exhibit an increase in the observed period that is more than 
the expected error with respect to Drake's mid-eclipse 
times. More observations are necessary to confirm or 
discard these trends. 

CSS07125, CSS080408, CSS09704, CSS09797, 
CSS21357, and CSS25601 wer e also disco v ered to be 
eclipsing WD+dM binaries by IDrake et all (|2010l ) with 
one mid-eclipse time in 2005 April or May listed, except in 
the case of CSS21357 where no eclipse time is presented. 
No further studies are available for these systems. 

PTFEBII.44I, PTFEB28.235 and PTFEB 28.852 are 
three eclipsing systems recently discover ed bv | Law et all 
(|2012f ) during the PTF/M-dwarf survey (|Law et al.ll2011l) 
for transiting planets around M-dwarfs. Only one mid- 
cclipse time is available for each system in 2010 August, 
September, and November, respectively. 

KIC-10544976 is a WD+MS PCEB in the field of 
Keple r mission, recently published by lAlmenara et al.l 
(|2012l ). The authors list then mid-eclipse times: six in 2005, 
three in 2006, and one in 2008. Currently, the system has 
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been continuously monitored by Kepler over two years, 
which makes it a good candidate to look for changes in 
the period in the near future. 

SDSS J0821+4559, SDSS J 0927+3329, SDSS 
30946+2030, SDSS J0957+3001, SDSS Jl 021+ 1744, 
SDSS J102S+0931, SDSS J1057+1307, SDSS J1223-0056, 
SDSS J 1307+2156, SDSS J 1408+2950, SDSS J HI 1+1028 
and SDSS J2235+U28 are all WD+MS PCEBs that 
were fo und to be eclipsing very recently bv iParsons et al.l 
2012al) bv correlating the SDSS WD+MS catalog 



Rebassa-Manserg as et al.l l2012f) with Catalina Real-time 



Transient Survey light curves. Only one mid-eclipse time 
is available for each system. 



A3. CVs with suspected planets 



an eclipsing polar (or 
by iBeuermann et al.l (fT98 



UZFor is 
discovered 

riod has been analyzed by several authors (e.g.. iRamsav 



AM Her type CV) 
The orbital pe- 



19941 H mamura fc Steiman-Camero n 1998; Perrvman et al 



2001), w ho detecte d varia tions from linearity in the O-C 
diagram. iDai et al.l (|2010ft detected an increase in the or- 
bital period and possible cyclical changes, which they sug- 
gest may be pro duced by a th i rd low -mass object in the 
system. Recently, iPotter et al.l (|2011l ) have combined new 
high-speed photometry spanning ten years with previous 
mid-eclipse times to obtain a baseline of 27 years. They 
find that the O-C diagram is described best by a combi- 
nation of two cyclic elliptical terms, probably due to the 
presence of two giant planets, and a secular variation. 
HUAqr is also an eclipsing polar discovered by 



Schwope et all (119931) during the ROSAT All Sky Survey. 
Schwarz et al.l (|2009h updated the ephemeris and listed 72 
eclipse egress times obtained between 1993 and 2007, re- 
vealing c omplex dev i ations from a linear trend in an O-C 
diagram. Qia n et all (|201lf ) combined these points with ten 
new eclipse egress times obtained between 2009 May and 
2010 May, and find two cyclic variations in the O-C curve 
and a long-term period decrease that cannot be explained 
by gravitational radiation or Applegate's mechanism. They 
propose that the cyclic variations observed are due to the 
presence of two giant planets, and the long-term period de- 
crease may reveal a third planet. However, the proposed 
orbital param eters for the two po ssible planets have been 
questioned by iHorner et al.l (|201ll) , who find the solutions 
extrem ely unstable on short time scales. IWittenmver et all 
(|2011l) used the same 82 eclipse egress times to fit different 
single- and double-planet models to the O-C diagram. They 
found that the best fits are obtained with two planets and 
that there is no need to invoke a third planet. Although 
the ir solutions are sig nificantly different from the one given 
by iQian et all (|201lh . the new parameters for the possi- 
ble planets are still dynamically unstable, which casts some 
doubt on their existence. They speculate about other mech- 
anisms being responsible for the observed variations, such 
as changes in the shape of the secondary star due to dynamo 
effect. However, in a very recent work, Go zdziewski et aLl 
(|2012ft publish almost 60 new eclipse egress times with bet- 
ter accuracy. Combining these with reanalyzed previous 
data, the authors find that a single circumbinary compan- 
ion gives the best explanation for the O-C curve. 

DP Leo was the first eclipsing polar discovered 
(|Biermann et al.lll985l). A decrea se in the binary per i od wa s 
noticed bv lSchwooe et all (|2002h and iPandel et all (|2002h . 



IQian et ail (|2010aft find a reversal of this trend, suggesting 
a sinusoidal varia tion that may be related to the presence 
of a giant planet. IBeuermann et all ([201 ll ) obtained accu- 
rate mid-eclipse times of the WD between 2009 March and 
2010 February, updated the ephemeris, and combined their 
results with all the mid-eclips e time s available by 2002, 
as published bv lSchwope et al.l (|2002lb The data spanning 
more than 30 years since 1979 suggest there is a third body 
orbiting the binary, most likely a giant planet. 
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